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Several microscopic pathways have been proposed to explain the large magnetic effects observed in organic 
semiconductors, but identifying and characterising which microscopic process actually influences the overall 
magnetic held response is challenging. Pulsed electrically-detected magnetic resonance provides an ideal plat¬ 
form for this task as it intrinsically monitors the charge carriers of interest and provides dynamical information 
which is inaccessible through conventional magnetoconductance measurements. Here we develop a general 
time domain theory to describe the spin-dependent reaction of exciton-charge complexes following the coher¬ 
ent manipulation of paramagnetic centers through electron spin resonance. A general Hamiltonian is treated, 
and it is shown that the transition frequencies and resonance positions of the exciton-polaron complex can be 
used to estimate inter-species coupling. This work also provides a general formalism for analysing multi-pulse 
experiments which can be used to extract relaxation and transport rates. 


I. INTRODUCTION 

Organic devices with non-magnetic electrodes may display 
large magnetoresistances at room temperaturJ^^. Magnetic 
fields of 10 mT or less are sufficient to induce large changes 
in sample conductivit}PI, with reports of changes exceeding 
2000 percent in one-dimensional wireP. There is interest in 
harnessing these field-induced effects to enhance technolo¬ 
gies including solar cellP, organic light-emitting diodeP, and 
cheap magnetic sensorP. This requires sophisticated meth¬ 
ods for characterizing and engineering of the underlying spin 
processes. 

Several detailed modelP® have been developed to explain 
the ubiquitous magnetoconductance effects observed in or¬ 
ganic devices. In recent years there has been significant 
progress in distinguishing potential spin processes through 
lineshape and voltage fingerprint^!^, however, these types of 
quasi-static measurements and models are fundamentally lim¬ 
ited: Dynamical rates and energetic couplings are obscured by 
the effects of spin ensemble averaging and the intrinsically in¬ 
coherent nature of the measurement. Coherent spectroscopic 
techniques, such as pulsed electron paramagnetic resonance 
(EPR), provide a way to independently identify and quan¬ 
tify the underlying physical processes because they allow sub¬ 
ensemble of spins to be selectively addressed and manipulated 
within the coherence 

Pulsed Electrically Detected Magnetic ResonancJ!^ 
(pEDMR) in particular provides a strong spectroscopic probe 
with which to investigate these proce sses, as different spin 
processes can readily be disentangled^^^^ and their impact 
on device operation directly observed from the transient 
dynamics. Pulsed methods are able to access coherent time 
domain information and by adjusting the applied pulse 
sequence different aspects of the spin ensemble can be 
probed, allowing a holistic picture of the spin physics to be 
developed and rigorously testecP^. Past studies have primarily 
focused on pristine Poly[2-methoxy-5-(2-ethylhexyloxy)-1,4- 
phenylenevinylene] (MEH-PPV) devices in which polaron 
pairs are the dominant spin-reaction pathway, and a detailed 
picture of the relevant energy and time scales in these devices 


has been developecJEC^. Dissociation, recombinat ion, and 
intersystem crossing rates have been characterizecP^^, as 
well as the strength of the electro-nuclear coupling and 
spin-spin interactions within the paii^S. Recent improve¬ 
ments in pEDMR measurement methods have substantially 
improved signal-to-nois^, and a variety of standard one and 
two-dimensional experiments from pEPP[!^ are now tech¬ 
nically feasible. In contrast, energy and time scales cannot 
be effectively separated in quasi-static magnetoconductance 
measurements; Ultra-small field effects are consistent with ei¬ 
ther a weak exchange interactiorP!! or slow carrier hoppin 
Differentiation betw een th ese two scenarios is possible with 
pulsed measurement however no comparable fingerprint 
is known for the static case. 

An extensive body of literature has been developed for de¬ 
scribing the time domain response of polaron pairs in mag¬ 
netic resonancJ^SEH^ while comparatively little is known in 
relation to possible three-particle processes. These are ex¬ 
pected to include triplet-exciton polaron quenching or the dy¬ 
namics of a bipolaron stabilized by a counter-ion in organic 
systemJ^. Other three-particle processes also include spin- 
dependent Auger recombination in nanocrystalJ^, trion states 
in silicorP^, and electrons interacting with a common nuclear 
batli^. 

In this work we develop a general formalism for describ¬ 
ing the spin-dependent response of three-particle complexes 
in pulsed magnetic resonance, thus providing the theoretical 
foundations for the quantitative characterization of the time 
and energy scales present in these systems. 

There are a number of different mechanisms invoked in 
the literature to describe the triplet-exciton polaron (TEP) 
quenching process. Eor instance, Desai et aP described a 
mechanism in which free polarons scatter from trapped triplet- 
excitons, causing the exciton to cross form the triplet to singlet 
manifold. 

Tl + D±i/2 ^ (Ti...D±i/2) -> D±i/2 + Sq (1) 

Where Ti is an excited triplet exciton, D is a free charge 
carrier and S is a singlet exciton. The reaction rate was as- 
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sumed to be spin-dependent, although no microscopic mech¬ 
anism was provided. In contrast, Koopmans and co-workers 
describe a quenching process which depends on the overall 
quartet or doublet content of the pair and may increase or 
decrease the polaron’s mobilit}^^. Other implementations in¬ 
clude a site-blocking mechanisnPSIand an exciton-polaron ex¬ 
change procesJ^. 

Similarly, there is substantial variation in the types of in¬ 
teractions assumed within the three-spin complex and the 
conditions under which they will be observed. The triplet- 
exciton polaron (TEP) process has been discussed in terms 
of weakly-interacting precursor pairJ^, strongly-bound trion 
stateJ^, and even bipolaron-counter ion complexeJ^. Each 
of these scenarios will produce a unique magnetic field- 
dependence and differentiation is an important step towards 
purposeful tuning of spintronic devices. The conditions un¬ 
der which the TEP mechanism will appear are also debated: 
High magnetic-field features from magnetoconductance mea¬ 
surements indicate that the triplet process is important at room 
temperaturJEl^ in contradiction with the sh^ decrease in 
triplet lifetimes observed at high temperature^ and in earlier 
pEDMR experiments^. 

In this paper we develop a time domain theory of an inter¬ 
acting spin-1-Hspin-1/2 pair which can react in proportion to 
its overall doublet (S=l/2) content, analogous to the theory 
of Ern and MerrifielcP^. The transition frequencies and reso¬ 
nance positions of an exciton-polaron complex are extracted 
from the canonical Hamiltonian, providing an analytic de¬ 
scription of the transition between a weakly spin-interacting 
exciton-polaron complex and a strongly-bound trion state. 
This result can be used to analyze pEDMR experiments on 
organic devices, allowing the intra-pair spin-spm coupling to 
be extracted, as demonstrated in related wor]p3. In this pa¬ 
per we concentrate on transitions following a single pulse, 
however the formalism and can be applied niore generally to 
multi-pulse sequences such as the Hahn echcPEll. As such, it 
provides the foundations for quantitative investigation of the 
triplet-exciton polaron quenching process with pEDMR. 


II. GENERAL FEATURES 

Long-lived triplet excitons can interact with trapped po- 
larons to form complexes, which may undergo a charge- 
reaction process which releases the trapped polaron and al¬ 
ters the free charge carrier density. Although a similar pro¬ 
cess is possible with singlet excitons, their density will be 
much loweil^ and we therefore ignore their contribution in 
this work. The spin-independent formation of exciton-charge 
complexes is likely to dominated by the reorganization energy 
associated with the shared lattice distortion^ due to the weak 
spin-interaction between the exciton and polaron. 

Eollowing its formation, an exciton-charge complex may 
i) dissociate back into its constituents or ii) recombine if it 
is in the doublet (spin-1/2) manifold, which will release the 
trapped polaron into the valence or conduction band, increas¬ 
ing the charge carrier mobility (see Eig. &)• The exact states 
involved in the reaction are unclear, however the weak spin¬ 


coupling seen experimentall}® suggests that a strongly-bound 
trion state is formed prior to recombination. 

[Exciton + Polaron] Exciton + Polaron (2) 

[Exciton + Polaron] [Trion] Sq Polar on"" (3) 

Spin resonance causes a mixing of the doublet and quartet 
manifolds, causing a net change in the recombination rate and 
average charge carrier mobility (conductivity). As in the two 
spin polaron-pair casJ^the second order contributions due to 
change in the electron and hole bath concentrations and mo¬ 
bilities are negligible, thus the expected changes in mobility 
are given by 


AlJ,e(t) = TL{Re{t) - Re,eq) 

AHhit) = TL{Rh{t) - Rh,eq) 

( 4 ) 

Where tl is the carrier lifetime, ne^h are the number of free 
electrons and hole, and R{t) and Req are the time-dependent 
and equilibrium recombination rates. 

The increase in charge carrier mobility from the charge re¬ 
action may yield either a current enhancement or quenching, 
depending on whether the majority or minority carriers are 
released. This counterintuitive behaviour arises due to the in¬ 
terplay between coulombic attraction and carrier transporP^. 
In general, we expect the charge reaction to release the major¬ 
ity carrier and provide a current enhancement. In the unipolar 
regime, which we treat as electron-rich, the change in current 
will be proportional to the number of released carrier. 

Acr(t) = eApe{t)ne (5) 

Eor a general description of the spin dynamics, we need to 
include the interaction of the spin ensemble with the environ¬ 
ment. This includes charge carrier generation, dissociation, 
recombination and relaxation. These infiuences may be in¬ 
cluded through the stochastic Liouville equatiorPI 

^=^-[p,H] + S[p] + R{p-po} (6) 

Where p is the density matrix and represents the occupa¬ 
tion probability of each state, H is the coherent spin Hamil¬ 
tonian and describes the exciton-polaron spin complex, S is a 
stochastic operator and represents the generation and loss of 
charge carriers, and R is the Redfield relaxation matrix and 
incorporates incoherent effects such as interactions with the 
surrounding phonon bath. 

A. Hamiltonian 

The spin Hamiltonian describes the coupling within the 
exciton-polaron complex and its interactions with the sur¬ 
rounding nuclear environment. As it completely encapsu¬ 
lates the coherent spin dynamics it can be used to derive the 
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FIG. 1. a) The general picture of the triplet exciton-polaron quenching mechanism. Exciton-polaron complexes are formed spin-independently 
due to the reorganization energy associated with a shared nuclear distortion. A complex may spin-independently dissociate back in a separate 
exciton and polaron or may undergo a exciton-charge reaction, possibly due to weak exciton-polaron exchange. Intersystem crossing will mix 
the doublet and quartet states prior to recombination or dissociation, b) A diagramatic representation of the Hamiltonian for an exciton-polaron 
complex with arbitrary coupling. The exciton will have internal exchange and dipolar coupling, and there may also be exchange and dipolar 
coupling between the exciton and the polaron. 


resonance positions and transition frequencies which allow 
the mechanism and coupling strengths to be identified. Both 
the exciton and polaron couple to the nearby nuclei through 
the hyperfine interaction and to the external magnetic field 
through Zeeman splitting. There is also dipolar coupling 
within the exciton, and dipolar and exchange coupling be¬ 
tween the exciton and polaron (see Fig. The exciton 

has a large internal exchange (J ^ 0.7 eV) which separates the 
singlet and triplet levels. As the internal exchange is signif¬ 
icantly larger than all other energy scales in a typical exper¬ 
imental regime, the singlet level becomes decoupled and as 
such the exchange interaction can be excluded from the effec¬ 
tive Hamiltonian. 


transform 

Henergy = HU (8) 

where 
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with 


H = PeB • qt ' St PdB ' qd ' Sd^ 

" -V-" -V-" 

Exciton Zeeman Polaron Zeeman 

E ln • An • /tz + / , ^Di • Adi • loi + 

v* v* 

* Exciton Hyperfine * Polaron Hyperfine 

St ' Dt ’ St “h JtdSt * Sjj 

■ v ■ v 

Exciton Dipolar Exciton-Polaron Exchange 

+ St • Dtd • Sd 

Exciton-Polaron Dipolar 

(7) 

An expanded version of this matrix is given in Appendix 
[A[ The commonly used high field approximations are made, 
which include an isotropic g-factor, isotropic hyperfine, and 
only retaining the secular terms of the dipolar couplinj^. 

Starting from the triplet-doublet (product) basis we can ro¬ 
tate the Hamiltonian to the energy eigenbasis using the unitary 


cot 20 = 


Dt 


^ + ^ - UJp+UJT 

+ V2Jtd 


and 


cot 2(j) = 


Dt 


^ ^ + wr, - wt 

+ V2Jtd 


( 10 ) 


( 11 ) 


This unitary transform preserves the physics of the Hamil¬ 
tonian, while simplifying analytic and numerical calculations. 
For instance, the transition frequencies simply correspond to 
the off-diagonal matrix elements. Importantly, the angles 0 
and p represent the degree of coupling between the exciton 
and polaron. Setting 0^(j) = {) corresponds to two independent 

particles and = cos“^ (the maximum rotation an¬ 
gle possible) corresponds to the formation of a spin-| trion. 
Due to the convenient measure of spin-spin coupling that 0 
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and 0 provide we will refer to them throughout the text as the 
spin mixing angles. 

The three natural basis sets for describing an exciton- 
polaron complex are the triplet-doublet (product) basis, the 
quartet-doublet basis and the energy eigenbasis. The transfor¬ 
mation between the triplet-doublet and quartet-doublet bases 
is given by the Clebsch-Gordan coefficients. 


Urot{t) = (18) 

In the rotating eigenbasis we have the Hq energy levels on 
the diagonal entries and the perturbative spin mixing Hi terms 
on the off-diagonal elements. 


\T+ l) = Q 



It is important to note that the term ‘doublet’ refers to a state 
with S=l/2, which may either be a single spin-1/2 particle or 
the I S', mg) = 1^, dz^) trion states. The transformation from 
the triplet-doublet to the energy eigenbasis is given by 
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Eq _ 


The dressed rotating basis states will have the same form as 
the bare energy eigenbasis states, and for most situations they 
may be treated equivalent!}!^. 

B. Limiting behaviour of the Hamiltonian 

While numerical calculations provide a general method 
for calculating the exciton-polaron dynamics from the spin 
Hamiltonian, this approach limit the physics insight which we 
can obtain. To understand the general properties of the system 
we now look at a number of limiting cases. 


Penergy p (12) 

with the triplet-doublet states given by 
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cos</>|r_ t)-cos(/)|ro 4 ,) 
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It is conceptually and computationally expedient to split the 
Hamiltonian into a large static component (Hq) comprising of 
spin-spin coupling and the static Zeeman splitting and a small, 
time-dependent driving term (Hi) which corresponds to the 
time-dependent Zeeman splitting. 


1. Weak coupling: 0,(j) ^ 0 

When the exciton-polaron coupling is extremely weak 
(6>, 0 ^0) the spins will be driven independently, and the 
eigenbasis is given trivially by the product states. The Hamil¬ 
tonian for this scenario is 
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H = Ho ^ Hi (15) 


Hi = 2Bi(e-^^^ -h + S^) (16) 

At high field the counter-propagating component will pro¬ 
duce no first-order effects, so we can move from the energy 
eigenbasis to the rotating frame and perform the rotating wave 
approximation. 

BH 

Hrot = UltHUrot - Ul^ — Urot (17) 


F+ti 
Tot 
T t 

Tot 

Lt- tJ 

The polaron resonances occur at a frequency of ^ 
with a Rabi frequency of ^Bi and the exciton will have res¬ 
onate at frequencies of + cct ± ^ with nutation frequen¬ 
cies of s/2^Bi. This is shown graphically in Fig. [^). There 
are three polaron transitions which occur independently of the 
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FIG. 2. The behaviour of an exciton-polaron complex depends strongly on the strength of the exciton-polaron coupling. Solid lines represent 
strongly-driven transitions and dashed lines represent transitions which become forbidden in the weak or strong coupling limits, a) In the weak 
regime the basis states are close to the product states and the transitions are the polaron transition, with a Rabi frequency of 7 B 1 , and 
the T+ o To, To o T_ exciton transitions with a Rabi frequency of V2jBi. The I T+ T- transition becomes weakly-allowed in the 
presence of non-zero coupling, b) In the intermediate coupling regimes the basis states are an admixture of the product and quartet-doublet 
states and up to eight transitions are visible with Rabi frequencies between 0 and VS'yBi . c) In the strongly-coupled regime quartet and 
doublets are formed. Strong driving occurs within the doublet and quartet manifolds and driving between the manifolds becomes arbitrarily 
slow, leading to a vanishing signal. The Rabi frequencies of these transitions are \/3^Bi and 2^Bi . 


exciton state, and likewise four possible exciton transitions 
which are independent of the polaron, shown by the solid 
lines. The single dashed line represents an Am^ = 1 tran¬ 
sition between the |t ^-) and T+) states which becomes 
weakly-allowed in the presence of non-zero exciton-polaron 
coupling. Eventually the spin-spin coupling becomes large 
compared to the hyperfine interaction and the transitions can 
no longer be described as excitonic or polaronic transitions, 
depicted in Fig. [^. 


2. 


Strong coupling: O^cf) cos 



When the exciton-polaron coupling is extremely strong 
(6>,0 ^ cos“^ ^ formed in the 

quartet-doublet basis. The exciton and polaron can no longer 
be considered separate entities and they will always nutate to¬ 
gether. The Hamiltonian is given by 


H = 
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( 22 ) 
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The mixing will be between the rris = ±3/2 and rris = 
±1/2 quartet manifolds with nutation frequencies of 
and mixing within the = ±1/2 doublet and quar¬ 
tet manifolds with Rabi frequencies of 2^Bi, as shown in 
Fig. 1^:. The four solid lines represent the strongly-allowed 
transitions of the trion state and the dashed lines represent 
quartet-doublet transitions which vanish in the strong cou¬ 
pling regime. Counter-intuitively the rris = ±3/2 quartet 
states actually couple to the rris = ±1/2 doublet states, which 
is nominally an = 2 transition. Analogous to the exciton 
half-field transition, it becomes weakly-allowed due to off- 
diagonal mixing in the quartet-doublet basis. Because the 
strongly-allowed transitions induce no overall change in the 
doublet or quartet content of the complex the trion states will 
be electrically undetectable. 

In this paper we analyze the transition between the uncou¬ 
pled and trionic regimes. 


C. Complex relaxation and recombination 

The recombination and dissociation rates of each eigenstate 
will be the weighted average of its quartet and doublet content. 
In analogy with the polaron paii^case we can form new rates 
for the admixture states. 


with the corresponding eigenstates 


n=TD\{i\D)\^ +rQ\m? 


(24) 
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di = dD\{i\D)\^ + dQ\{i\Q)f 


(25) 


occupation of each state corresponds to when the pair genera¬ 
tion and loss rates are equal. 


The spin loss rate is the sum of the dissociation and recom¬ 
bination rates. 


+ di 


(26) 


Pi = 


9i 

Vi + di 


(30) 


When doublet states are lost much faster than quartet states 
there will be a strong polarization of pure quartet states. 


A spin-dependent depopulation of the doublet-state has 
been observed in other systems and is consistent with sim¬ 
ple energetic arguments and spin-selection ruleJ^. In the 
weak coupling limit {0,(j) ^0) the T_ i states have 

I doublet-content, and the Tq states each have ^ doublet- 
content. 


\rD H- \rQ 

Itd + Irq 

fD + |rg 
^VD + srg 

L rg J 


(27) 


In the opposite, strong coupling limit (O^cf) ^ cos 
the states become pure quartet and doublet states. 



td 

td 

SQ. 


(28) 


Between these two limits there will be a monotonic change 
in the recombination rates as the product basis morphs into the 
quartet-doublet basis. 


P0 = i(T+t+T_;) (31) 

In organic devices the polarization effects arise primar¬ 
ily from the different rates of charge carrier losJI^. In con¬ 
trast, porphyrin system may generate a polarization within the 
triplet-doublet pairs from the intersystem crossing of light- 
induced singlet-doublet stateJ^^^il 


E. Transient readout 

The real-time detection of sample conductivity changes 
during resonant excitation is technologically challenging^. 
Consequently, experiments are usually performed using a 
pump-probe approach, which involves monitoring the (slow) 
response of the system following the (rapid) magnetic reso¬ 
nant excitatioiP^. 

In the triplet-exciton polaron quenching mechanism recom¬ 
bination will trap or detrap the polaron and change the average 
charge carrier mobility. This manifests macroscopically as a 
change in the sample conductivity. We can therefore relate 
the recombination rate or overall doublet content of the spin 
ensemble in our model to an experimental change in conduc¬ 
tivity. Larmor beating will dephase the exciton and polaron 
on a time scale much faster than the transient readout, and we 
can therefore treat the spin ensemble incoherently using rate 
equationJ^ 


D. Complex formation and polarization 


(32) 

i 


We can’t discuss the transitions of an exciton-polaron com¬ 
plex without discussing their formation and steady state po¬ 
larization: The presence of an electrically-detectable signal 
intrinsically requires a net transfer of the spin population. 

Exciton-polaron complexes are believed to form spin- 
independently. A spin-dependent formation has been pro¬ 
posed in the context of polaron pairs, and a similar effect for 
exciton-polarons will change the intensity of the transitions, 
but leave the qualitative conclusions of this study unaffected. 
Regardless, the generation of each state needs to be weighted 
by a Boltzmann factor to account for the thermal distribution. 


9i= 9 


exp(-i^) 


(29) 


At high temperatures polarization effects can be ignored 
and a uniform generation rate can be used. The steady state 


Where Spa is the change in occupation probability and the 
recombination rate of each state. The overall change in cur¬ 
rent induced by each transition will depend on the precise 
choice of integration window and dissociation, recombina¬ 
tion, and intersystem crossing rates. In this work, we are con¬ 
cerned with the resonance positions and transition frequen¬ 
cies, which are unaffected by these incoherent rates. Nonethe¬ 
less, we sketch how quantitative predictions may be made 
by discussing the minimalist scenario in which intersystem 
crossing and recombination are slow compared to dissocia¬ 
tion, which produces an observable which is proportional to 
the doublet content of the occupied states. We also assume 
that all recombination events are recorded. 

Q = y] / riexp (-(rj + di)t)dt ~ E] ^ ~ X] I {D\(i>i) P 

i Jo , «i 
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Where is the i-th eigenvector and D is the doublet man¬ 
ifold. An unknown scaling factor is required to relate Q, the 
overall microscopic change in charge carrier mobility, to a real 
change in sample conductivity. 


III. TRANSITION ANALYSIS 

In the standard implementation of a pEDMR experiment 
the transition frequencies are mapped as a function of energy 
by sweeping the magnetic field with a constant excitation fre¬ 
quency and measuring the change in charge as a function of 
the microwave pulse length. This may result in a number 
of different resonances each with a particular transition fre¬ 
quency. These serve as a ‘Rabi fingerprint’ which can be used 
to identify and quantify the types of coupling present. We now 
calculate these ‘fingerprints’ directly from the spin Hamilto¬ 
nian. 


A. Resonance position 


The photon energy required to resonantly drive a transi¬ 
tion is equal to the energy difference between the two eigen¬ 
states of interest. The energies of the Hamiltonian described 


by equation (A7) which are used to calculate the resonance 
positions are given in Table 

In the weak coupling limit the three polaron transitions will 
occur at the Larmor precession frequency and the four exciton 
transitions will be paired and will be offset the intra-exciton 
dipolar coupling. The multiplicity of each transition is due 
to possible states of its partner, so the exciton has identical 
transitions whether the polaron is in the |t) or ||) state. In the 
strong-coupling the rris = 1/2 transitions {Qi ^ ^ 

D_i) will occur at the Larmor precession frequency and the 
other quartet resonances (Q±3 ^ ^ large 

separation from the substantial zero-field splitting. 


B. Transition frequencies 

The frequency of the Rabi oscillations are indicative of the 
types of states and coupling that are present, and are calcu¬ 
lated by solving the 2 x 2 eigenvalue problem of the relevant 
subspace. The frequencies of the Rabi oscillations have the 
form 


Weak coupling 


Strong coupling 



FIG. 3. The transition frequencies of an exciton-polaron complex 
as a function of the mixing angle. In the uncoupled regime {0 = 0) 
the polaron will nutate at a frequency of 7 B 1 and the exciton will 
nutate at a frequency of . Weak coupling splits the exciton 

and polaron resonances and permits the f T+ of T_ transition to be 
weakly driven. In the strong coupling limit pure quartet and doublet 
states are formed. Quartet mixing will occur at frequencies of 
(Qi/2 ^ Q-1/2) and \/37Bi (Q±3/2 ^ Q±i/2)- Driving can also 
occur within the doublet manifold at a frequency of jBi . Driving 
between the doublet and quartet states becomes arbitrarily slow in 
the strong coupling limit. Transition frequencies are calculated by 
using equation 


The transient nutation frequency, aij, for each possible 
transition is given in Table |T| and are shown graphically in Fig. 

as a function of the mixing angle 0 (with 0 = (/)). It can 
be seen that the nutation frequencies change distinctively as a 
function of mixing angle, and this can be used to differenti¬ 
ate triplet-exciton polaron quenching from other mechanisms 
such as the triplet-triplet annihilation, as well as identifying 
the relevant coupling regime. 


TABLE 1. The transitions frequencies of an exciton-polaron complex 
as a function of the spin mixing angles 


Transition 

Frequency ( 7 B 1 ) 

(l)-( 2 ) 

y/2 cos 0 + sin 0 

(l)-(4) 

cos 0 — \/2 sin 0 

(2)-(3) 

a/2 cos 0 cos 4 > — sin 0(cos 0 + a/ 2 sin 0) 

(2)-(5) 

cos (/)(cos 0 + a/ 2 sin 0) + a/2 cos 0 sin 0 

(3)-(4) 

sin ^(sin cj) — a/2 cos 0) — a/2 cos 0 sin 0 

(3)-(6) 

cos (j) — a/ 2 sin cj) 

(4)-(5) 

— cos 0(sin 0 — a/2 cos 0) — a/ 2 sin 0 sin 0 

(5)-(6) 

sin (j) + a/2 cos cj) 


flij = + (w - (34) 

Where aij is the frequency coefficient which may be read 
from the diagonal matrix entries, uo is the applied photon 
frequency, and ujij is the resonant (i,j) transition frequency. 
When a resonant excitation (cc = uJij ) is applied the nutation 
frequency reduces to 

^ij = aijjBi (35) 


C. Regimes 

We are able to identify several driving regimes depending 
on the relative energy scales present. Figure shows the tran¬ 
sition frequencies as a function of mixing angle. Additional 
regimes are possible when the exchange and dipolar coupling 
are able to cancel out on either the diagonal or off-diagonal 
elements of the Hamiltonian. 
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Uncoupled regime: The simplest case, in which there 
is no interaction between the exciton and polaron 
(Jtd, Dtd = 0)- The polaron will nutate at jBi and 
the exciton will be split into two lines with transition 
frequencies of When the intra-exciton dipo¬ 

lar vanishes the exciton levels become harmonically 
spaced and a multi-transition analysis is required, which 
we do not treat in this paper. 

Effectively uncoupled regime: Similar behaviour occurs in 
the effectively uncoupled regime, which occurs when 
\Dtd H- 3Jtd| << Dt. The off-diagonal mixing van¬ 
ishes and the physics will be identical to the weakly- 
coupled regime except that the transitions have a greater 
splitting. 

Weakly coupled regime: The generic weakly coupled 
regime corresponds to when the exciton-polaron 
coupling is weak {\Dtd\i \Jtd\ « Dt). The exciton 
and polaron resonances become weakly correlated - the 
and To(t^i) transitions will have slightly 
different transition frequencies, for instance, and the 
T+ T_ t transition becomes weakly allowed and 
will provide a low frequency background. 

Isoenergtic regime: An interesting case occurs when 
Dtd = \Jtd and Jtd « Dt- The energy levels 
will be unshifted to first order, but there will be a small 
shift in the transition frequencies and the effective 
recombination rates. 

Trionic regime: The trionic regime will occur when there 
is strong coupling between the exciton and polaron 
states from exchange or dipolar coupling. In either 
case, the exciton and polaron become strongly entan¬ 
gled and form either quartet or doublet states. Driving 
between the quartet and doublet states (Qj. s ^ i, 
Q^i ^ D^i) will have a vanishing small probabil¬ 
ity. Driving will also occur between the doublet states 
(Di ^ D_0 at a frequency of jBi and between the 
quartet states (Qj-s ^ ^ 

quencies of and 2^Bi, respectively. As these 

transitions leave the overall quartet-doublet content un¬ 
changed, they will not produce an observable signal in 
the absence of spin-orbit coupling or some other addi¬ 
tional spin-reaction pathway. 

Intermediate coupling regime: Between these two extremes 
lies the intermediate coupling regime which involves 
neither pure triplet-doublet, nor quartet-doublet basis 
states. There will eight possible transitions with fre¬ 
quencies between zero and three ^Bi . Due to the large 
parameter space this regime encompasses a numerical 
approach is usually required. 

We have so far dealt with an exciton interacting with a po¬ 
laron, which corresponds to positive exchange and dipolar 
coupling using our definition. A bipolaron interacting with 
a counter-ion can be treated by including negative exchange 



FIG. 4. The transition frequencies of a bipolaron counter-ion com¬ 
plex, which is treated by introducing a negative mixing angle. At an 

angle of bipolaron is broken apart by 

the strong coulombic attraction from the counter-ion. If the coupling 
is increased further an exciton-polaron complex is formed for a mix¬ 
ing angle of ^2 = cos“^ cos“^ From mixing angles of ^2 

to O3 = cos“^ the transition frequencies of the exciton-polaron 
complex evolve in an identical manner to the positive mixing angle 
regime. Transition frequencies are calculated from equation ( [^ . 


and dipolar coupling, a scenario advocated by Shinai^^. Once 
again we can precisely extra the spin-spin coupling from the 
‘Rabi fingerprint’. Weaky coupling produces individually ad¬ 
dressable states with characteristic frequencies of jBi and 
\/2^Bi . As the coupling strength is increased the bipolaron 
will be tom apart and an exciton-polaron will be formed. Af¬ 
ter this point the states may form a trion, proceeding in a sim¬ 
ilar manner to above. The bipolaron can be considered to be 
separated for a mixing angle of 

— 1 r? —1 /l 

COS \ % — cos \ k 

e = - ^ (36) 

which produces transitions of ^and | 7 ^ 1 . At a mix¬ 
ing angle of 


6 » = cos ^ - cos ^ (37) 

an uncoupled exciton-polaron complex will be formed, 
with the expected transition frequencies of 7^1 and \/2^Bi 
for the polaron and exciton respectively. The infiuence of 
spin-spin coupling on the bipolaron is shown in Fig. 


IV. HALF FIELD 

It is possible to directly excite a Am^ = 2 spin-flip transi¬ 
tion between the T_ and T+ states of an exciton at the half 
field resonance. 

^ = 27 B 1 [ 57 , + 57 , + 57](e2-°* + e- 2 - 0 *) (38) 
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This transition is nominally forbidden, but becomes weakly 
allowed due to off-diagonal zero-field splitting which sup¬ 
presses spin as a good quantum number and creates admix¬ 
tures of the T+ and T_ states. 

To perform quantitative analysis of the half field transitions 
it is necessary to move to the doubly rotating frame and ex¬ 
pand the time-averaged Hamiltonian to second ordei^. 


[ ° ° [ \H{t2),H{h)]dhdt 

Jo Jo Jo 

This produces a term proportional to 


H = — 

Tn 


2 

(39) 


BiD sin 20 
^0 




(40) 


which is responsible for the Arris = 2 transition. To ob¬ 
tain the correct transition probability, it is necessary to inte¬ 
grate over all the possible orientations with the appropriate 
weighting. The ensemble half-field Rabi frequency will scale 
with the driving field and dipolar strength in the perturbative 
regime. 


V. MICROSCOPIC BASIS FOR RECOMBINATION 

So far we have discussed the spin-dependent recombination 
of exciton-polaron complexes. An obvious argument to ac¬ 
count for the different rates is spin-conservation: The ground 
state is S=l/2, so we should expect the doublet states to decay 
faster than the quartet states, however this ignores that both 
are orthogonal to the ground state: 


TABLE 11. The unnormalized eigenkets of a weakly coupled exciton- 
polaron complex when there is a small difference in the exchange 
coupling. The resultant mixing between singlet and triplet manifolds 
allows recombination of the doublet states to the singlet ground state. 


No Exchange 

Exchange 

tT+ 

tT+ 


4 5+ To - JJAJxd I 

6 "'6 

t^o 

4 To - ^^5 

\T- 

4 T- - 1 5- 

;r+ 

4 T+ + J’_ 

iS 

4 5+^f/-|,4To+^f/-|,4T+ 

4 To 

4 To + 4 5 

J -\- 3 

4T_ 

4T_ 


AJtD — Se^ex ' Se,pol Sh,ex ' Se,pol (44) 

The asymmetric exciton-polaron exchange interaction 
(AJtd) is isotropic and will be small in comparison with 
the exciton exchange (Jt). The new eigenstates may be ef¬ 
ficiently calculated analytically or numerically, but are gener¬ 
ally cumbersome to deal with. A fruitful path is to examine 
the role of exciton-polaron exchange through first order per¬ 
turbation theor})^ 


An = (no|AJTD|^o) (45) 


(Q|o,i) = o 


(41) 


|n) = + 






( 0 ) 


E 


( 0 ) 


(46) 


(Z)|0, i) = 0 (42) 

We propose a possible mechanism in which spin conser¬ 
vation holds, yet allows spin angular momentum to be trans¬ 
ferred between the exciton and polaron: A Dexter exchange- 
based pathway reliant on a small asymmetry in the exciton- 
polaron exchange. In the canonical Hamiltonian it is tacitly 
assumed that the exchange between the lone polaron and each 
of the constituent polarons are equal. Due to the opposite po¬ 
larity of the exciton polarons, one would expect there to be a 
small asymmetry in the two coupling strength^ 


Where is the energy shift and \n) are the new eigen¬ 
states in terms of the uncoupled states. The resulting eigenkets 
are given in Table [V| 

In this mechanism the exciton-polaron recombination oc¬ 
curs due to the doublet states becoming mixed with the singlet 
exciton. The rate of the reaction will depend on the recom¬ 
bination rate of the singlet exciton and the intra-exciton and 
exciton-polaron exchange strengths. 


{S\D) 


AJtd 

J 


(47) 


JTD,tot = Jtd + AJtd (43) 


R 


exciton-polaron 


= R 


singlet 


AJtd 

J 


(48) 


Where Jtd, tot is the total exchange coupling between the 
exciton and polaron, Jtd is the sum of the two microscopic 
exchange interaction, and AJtd is a small asymmetry be¬ 
tween the electron-polaron and hole-polaron exchange coal¬ 
ings. 


Since the singlet admixture will be proportional to the dou¬ 
blet content it provide the same dynamics as a process medi¬ 
ated directly via doublet content. A similar argument can be 
made involving recombination due to a slightly asymmetric 
dipolar coupling between the exciton and polaron. 
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VI. CONCLUSION 


Triplet-exciton polaron quenching plays an important role 
in the large magn etic field effects observed in organic 
semiconductori^ll^ We have developed a general time- 
domain theory which quantitatively describes the changes in 
samples conductivity due to this mechanism in pulsed Elec¬ 
trically Detected Magnetic Resonance. Population mixing 
due to resonant excitation causes a net change in the exciton- 
charge reaction rate, and hence the free charge carrier density. 

In particular, we have derived transition frequencies and 
resonance positions for an exciton-polaron complex with arbi¬ 
trary coupling. Our modeling indicates that there is a clearly 
discernible transition between an uncoupled exciton-polaron 


state and a strongly-bound trion. These situations may be dif¬ 
ferentiated through their transition frequencies and resonance 
positions. Our formalism provides a general basis for quan¬ 
titative analysis of the triplet-exciton polaron process through 
pEDMR. In the future multi-pulse schemes such as the Hahn 
echo and inversion recovery could be used to investigate re¬ 
laxation processes directly. 

VII. ACKNOWLEDGMENTS 

This work was supported by the Australian Research 
Council (ARC) (DP120102888 and CEl 10001027). TLK 
is supported by the Australian Renewable Energy Agency 
(ARENA)(7-S009). DRM acknowledges an ARC Euture Eel- 
lowship (ET130100214). 


* t.keevers@student.unsw.edu.au 

^ T. Francis, O. Mermer, G. Veeraraghavan, and M. Wohlgenannt, 
New Journal of Physics 6, 185 (2004). 

^ O. Mermer, G. Veeraraghavan, T. Francis, Y. Sheng, D. Nguyen, 
M. Wohlgenannt, A. Kohler, M. Al-Suti, and M. Khan, Physical 
Review B 72 , 205202 (2005). 

^ R. Mahato, H. Liilf, M. Siekman, S. Kersten, P. Bobbert, 
M. de Jong, L. De Cola, and W. van der Wiel, Science 341 , 257 
(2013). 

H. Hoppe and N. S. Sariciftci, Journal of Materials Research 19 , 
1924 (2004). 

^ A. Dodabalapur, Solid State Communications 102 , 259 (1997). 

^ W. Baker, K. Ambal, D. Waters, R. Baarda, H. Morishita, K. van 
Schooten, D. R. McCamey, J. M. Lupton, and C. Boehme, Nature 
Communications 3, 898 (2012). 

^ V. Prigodin, J. Bergeson, D. Lincoln, and A. Epstein, Synthetic 
Metals 156 , 757 (2006). 

^ P. Desai, P. Shakya, T. Kreouzis, and W. Gillin, Physical Review 
B 76 , 235202 (2007). 

^ P Bobbert, T. Nguyen, F. Van Oost, v. B. Koopmans, and 
M. Wohlgenannt, Physical Review Letters 99, 216801 (2007). 

P. Janssen, M. Cox, S. Wouters, M. Kemerink, M. Wienk, and 

B. Koopmans, Nature Communications 4, 2286 (2013). 

C. Boehme and J. M. Lupton, Nature Nanotechnology 8, 612 
(2013). 

C. Boehme and K. Lips, Physical Review B 24, 245105 (2003). 

T. Keevers, A. Danos, T. Schmidt, and D. R. McCamey, Nature 
Nanotechnology 8, 886 (2013). 

A. Schweiger and G. Jeschke, Principles of pulse electron para¬ 
magnetic resonance (Oxford University Press, 2001). 

D. R. McCamey, K. J. van Schooten, W. J. Baker, S.-Y. Lee, S.-Y. 
Paik, J. M. Lupton, and C. Boehme, Physical Review Letters 104 , 
017601 (2010). 

W. Baker, T. Keevers, J. Lupton, D. McCamey, and C. Boehme, 
Physical Review Letters 108 , 267601 (2012). 

D. McCamey, H. Seipel, S. Paik, M. Walter, N. Borys, J. Lupton, 
and C. Boehme, Nature Materials 7, 723 (2008). 

D. McCamey, S.-Y. Lee, S.-Y. Paik, J. Lupton, and C. Boehme, 
Physical Review B 82 , 125206 (2010). 

H. Malissa, M. Kavand, D. Waters, K. van Schooten, P Burn, 
Z. Vardeny, B. Saam, J. Lupton, and C. Boehme, Science 345, 
1487 (2014). 


F. Hoehne, L. Dreher, J. Behrends, M. Fehr, H. Huebl, K. Lips, 
A. Schnegg, M. Suckert, M. Stutzmann, and M. S. Brandt, Review 
of Scientific Instruments 83, 043907 (2012). 

T. Nguyen, B. Gautam, E. Ehrenfreund, and Z. V. Vardeny, Phys¬ 
ical Review Letters 105 , 166804 (2010). 

A. Schellekens, W. Wagemans, S. Kersten, P. Bobbert, and 

B. Koopmans, Physical Review B 84, 075204 (2011). 

V. Rajevac, C. Boehme, C. Michel, A. Gliesche, K. Lips, S. Bara- 
novskii, and P. Thomas, Physical Review B 74 , 245206 (2006). 

A. Gliesche, C. Michel, V. Rajevac, K. Lips, S. Baranovskii, 
F. Gebhard, and C. Boehme, Physical Review B 77, 245206 
(2008). 

M. Limes, J. Wang, W. Baker, S.-Y. Lee, B. Saam, and 

C. Boehme, Physical Review B 87 , 165204 (2013). 

J. Shinar, Laser & Photonics Reviews 6, 767 (2012). 

K. J. van Schooten, C. Boehme, and J. M. Lupton, 
ChemPhysChem (2014). 

A. Belolipetskiy, O. Gusev, A. Dmitriev, E. Terukov, and 

I. Yassievich, Semiconductors 48, 235 (2014). 

S. Weber, G. Kothe, and J. R. Norris, The Journal of chemical 
physics 106 , 6248 (1997). 

J. Song, N. Stingelin, A. Drew, T. Kreouzis, and W. Gillin, Physi¬ 
cal Review B 82 , 085205 (2010). 

M. Cox, P. Janssen, F. Zhu, and B. Koopmans, Physical Review B 
88 , 035202 (2013). 

M. Baldo and S. Forrest, Physical Review B 62, 10958 (2000). 

W. Baker, D. McCamey, K. van Schooten, J. Lupton, and 
C. Boehme, Physical Review B 84 , 165205 (2011). 

V. Ern and R. Merrifield, Physical Review Letters 21, 609 (1968). 

W. Baker, T. Keevers, C. Boehme, and D. R. McCamey, submitted 
to Physical Review Letters (2014). 

H. Huebl, F. Hoehne, B. Grolik, A. R. Stegner, M. Stutzmann, and 

M. S. Brandt, Phys. Rev. Lett. 100 , 177602 (2008). 

W. J. Baker, T. L. Keevers, J. M. Lupton, D. R. McCamey, and 
C. Boehme, Physical Review Letters 108 , 267601 (2012). 

S. StafStrom, Chemical Society Reviews 39, 2484 (2010). 

F. Bloom, M. Kemerink, W Wagemans, and B. Koopmans, Phys¬ 
ical Review Letters 103 , 066601 (2009). 

N. Atherton, Principles of electron paramagnetic resonance 
(1993). 

M. Asano-Someda, A. van der Est, U. Kruger, D. Stehlik, 
Y. Kaizu, and H. Levanon, The Journal of Physical Chemistry A 




11 


103 , 6704 (1999). 

Y. Kandrashkin and A. van der Est, Chemical Physics Letters 379 , 
574 (2003). 

A. van der Est, M. Asano-Someda, P. Ragogna, and Y. Kaizu, The 
Journal of Physical Chemistry A 106 , 8531 (2002). 

Y. Kobori, K. Takeda, K. Tsuji, A. Kawai, and K. Obi, The Journal 
of Physical Chemistry A 102 , 5160 (1998). 

F. Hoehne, C. Huck, M. S. Brandt, and H. Huebl, Physical Review 
B 89, 161305 (2014). 

J. J. Sakurai and S. F. Tuan, Modern quantum mechanics, vol. 104 
(Addison-Wesley Reading (Mass.), 1994). 

I 


L. Swanson, J. Shinar, and K. Yoshino, Physical Review Letters 
65 , 1140(1990). 


Appendix A: Hamiltonian 

The expanded exciton-polaron Hamiltonian defined is 
equation (|7]) is 


'El 

0 

0 

0 

0 

0 


OJp 

2 


\/2Dtd 

6 


2Dt 

3 


J 

x/2 


Dt 

2 ' 3 


3 

0 

V^Dtd 

6 

0 


_L _1_ 

' 2 ' 2 


■ (jJj^ 


J 


0 

_ ^/2Dtd J 

6 x/2 

0 

0 

0 


0 


0 
0 

A 0 

0 


UJs _ 2Dt _ U; £) Q 

' 2 3 2 ^ 

0 Eej 

(Al) 


r 


Where ^ 


— i + UJt ^ttd Eq — 


2 ' 3 ' 3 

_ 3cc>g I Dt I Dtp _ J_ _ ^p _ , , 

2 ' 3 ' 3 2 2 • 

The static Hamiltonian (Hq) may be transformed to the en¬ 
ergy eigenbasis using the unitary matrix from equation (|^. 


The perturbative Hi driving terms are strictly off-diagonal 
elements, which mix the eigenstates, and determine the tran¬ 
sition frequencies. 


3ujs . Dt Dtd j . . 

El — -^- - -^- - -— H- - -h iOT 


ujg J Dtd Dt ojt 


(A2) 


E2 = ^E- 


2 4 6 6 2 

{-{6 Dt — 2Dtd + 3J)^ + - {Dtd + 3J)^) (A3) 


_ ^ J _ Dtp _ Dt^ ^ ^ 

^ ~ ~2 ^ 6 

+ (l(6r>T - 2Dtd + iJ? + \{Dtd + 3J)2) (A4) 


^ ~ ~Y ^4 6 ^ ^ 

+ (-(fi^T — ‘^Dtd + 3J)^ + - {Dtd + 3J)^) (A5) 


_ _ Re I ee 

^ 2'4 6 6'2 

- {^{GDt - 2Dtd + iJ? + \{Dtd + iJf) (A6) 


Bi 2 = B 21 = (sin 0 -1- s/2 cos 0) 

Bi 4 = B/ii = (cos 0 — V 2 sin 0) 

B 23 = B 32 = —^(cos^sin^—\/2cos^cos0-h\/2sin^sin0) 
B 23 = B 32 = — (cos6>cos0+\/2cos6>sin0+\/2sin6>cos(/)) 
B34 = B43 = —^(v^cos^sin^—sin^sin 0 -hv^cos 0 sin^) 
-B36 = -Bes = (cose/) - V 2 sin^) 

B45 = B3/1 = —^(cos^cos0—v^cos^sin0+v^sin^sin0) 


-B56 = -Bes = ^(sin(/) + V 2 cos^!')) 


(A8) 


Removing the Zeeman contribution from the static mag¬ 
netic field (Bq) and substituting into equation provides 
a general description of the coherent dynamics of an exciton- 
polaron complex. 


Appendix B: Resonance Position 


3iUs . Dt Dtd J lod 


The resonance positions are given by the energy differences 
of the eigenstates. 
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TABLE III. The resonance positions for an exciton-polaron with ar¬ 
bitrary spin-spin coupling, calculated from equation ( |A7| ). In the ta¬ 
ble below X = ( I (6Dt — 2Dtd + + | {Dtd + 


Transition 

AE 

(l)-(2) 

JJt 3J I UJd I 1 V 

2 ' 2 ' 2 ' 2 ^ 

(l)-(4) 

Dt Djpd 3J I UJJJ 1 ujt V 

2 ' 2 4 ' 2 ' 2 ~ ^ 

(2)-(3) 

lot - 2X 

(2)-(5) 

lot 

(3)-(4) 

(jOt 

(3)-(6) 

3J Dtd Drjp . u)]j 1 COT 1 V 

4 2 2 ~'~ 2 ~'~ 2 '^^ 

(4)-(5) 

UJt + 2X 

(5)-(6) 

3J D'xj^ Dj- 1 ujqp 

4 2 2 2 ~ ^ 


■Ei 


Ej 


(Bl) 


These are calculated explicitly from equation ( |A7| ) and 
given in Table 

All of the transitions have a resonance position which de¬ 
pends on the spin coupling term ‘X’, which indicates that 
all of the transitions will become very broad as the exciton- 
polaron coupling becomes large. 


Appendix C: Microscopic basis for recombination 


The expanded (anti-symmetric) exchange matrix is 


"0 0 0 0 0 0 0 0 " 

0 0 0 0 0 0 
0%^00 0 0 00 
00000^00 
= Q 0 0 0 0 0 0 

0 0 0 ^0 0 0 

0 0 0 0 0 0 0 
00000 0 00 

(Cl) 

Where the coupling term is written in the expanded 
singlet/triplet-doublet basis. All of the matrix lie on off- 
diagonal elements, which means the first order effect will be 
to slightly admix the eigenstates, but not change the energies. 
This term therefore allows the canonical triplet-doublet states 
to project onto the singlet-doublet states. A similar argument 
can be made for the exciton-polaron dipolar coupling term. 










